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1. Introduction 



In a recent paper published in this journal , Di Francesco, Gaudin, Itzykson, and 
Lesage have studied a classical two-dimensional one-component plasma submitted to a 
quadrupolar electric field. The present paper is a supplement to the work of these authors. 

The model under consideration consists of N classical particles of charge e immersed 
in a uniform neutralizing background. The particles interact through the two-dimensional 
Coulomb potential — (e 2 /2n) In \z{ — Zj\, where the position of the i-th particle is defined 
by the complex number zf, there are also background particle and background-background 
interactions. A dimensionless coupling constant is (3 — e 2 /27r/cBT, where T is the temper- 
ature and ks is Boltzmann's constant. For the special value (3 = 2, the model is exactly 
solvable in a variety of geometries 2-5 . A simple case is when the system is a disk; then the 
background generates a one-body potential of circular symmetry. However, Di Francesco 
et al have considered the more general case when the one-body potential V(z, z) also has 
a quadrupolar component ; up to some additive constant, 



With this choice, the model can still be solved for (3 = 2. Di Francesco et al have studied the 
density of the A^-particle system in the large- N limit, and they found that the particles form 
a blob of constant density 1/tt, having the shape of an ellipse with semi-axes a = y/Ne' 1 
and b = y/Ne~^. 

In the present paper, we revisit the A^-particle system in a quadrupolar field. In Section 
2, we retrieve the elliptic shape of the particle blob by a simple macroscopic argument. 
In Section 3, we revisit the microscopic approach by using the analogy with a system of 
independent fermions in a magnetic field. In Section 4, we study the correlations near the 
boundary of the blob. 

2. Macroscopic equilibrium 





The one-body potential (1.1) can be thought of as generated by a uniform background 
with the shape of a large ellipse. Indeed, if the background charge density is — e/ix and the 
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ellipse has semi-axes A and B such that B/A = e~ 2At , a simple calculation 6 shows that, at 
a point z inside the ellipse, this background generates an electric potential V/e given by 
(1.1), up to some additive constant (here ksT = e 2 /4n). 

Let us now put N particles in, assuming N << AB (this condition will ensure that 
the particles stay well inside the background) . The particles will form a blob, and macro- 
scopic electrostatics requires that the total electrostatic potential be a constant inside that 
conducting blob. Therefore, the density and shape of the blob should be such that it 
generates an electric potential —V/e, up to an additive constant. The same calculation as 
above says that the blob which achieves that has the shape of an ellipse, with semi-axes 
a and b such that b/a = e _2M , and a constant particle density N/rcab = 1/ir. Thus, one 
retrieves the result of Di Francesco et al : a = y/We^, b = \TNe~ 11 . 



3. Magnetic picture 

The Boltzmann factor of the one-component plasma is, up to a multiplicative constant, 

N 

ex P { E I - z & + 2 ( tan M(4 + } l A ! 2 

where 



A = det 



1 Z\ Z-^ . . . Z^ 
1 Z 2 Zn . . . Z 2 



1 -v -v2 ~N-1 
1 Z]\f z N ... z N 



i>3 



Di Francesco et al have expressed this Boltzmann factor in terms of orthogonal functions. 
We now rephrase their work, in a way that we believe to be more transparent, using the 
magnetic language. Indeed, in the special case of circular symmetry (fx = 0), (3.1) was 
introduced as a squared wavefunction of a quantum system : iV independent fermions in 
a plane, interacting with a uniform magnetic field normal to the plane. Let us first review 
this picture 7 for \i = 0. 

The quantum Hamiltonian for a particle of mass 1 and charge 1 in a magnetic field is 
H = (l/2)(p — A) 2 , where p is the canonical momentum and A the vector potential. For 
describing a uniform magnetic field normal to the plane (the value of which is chosen as 



2), a possible gauge is A x = —y, A y = x (x and y are Cartesian coordinates in the plane; 
z = x + iy) . In terms of annihilation and creation operators for "left quanta" 7 

a = \i x + W + iPx ~Py) = ^{x + iy + d x + id y ) = + d- z 
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a + = -(x -iy- %Vx ~P y ) = ~(x - iy - d x + id y ) = -z - d z 

the Hamiltonian is a 

H = 2a+a+ 1. 

For obtaining a complete set of commuting observables, we must also introduce annihilation 
and creation operators for "right quanta" 

b= -(x - iy + ip x + p y ) = -(x - iy + d x - id y ) = ^z + d z 
b + = -(x + iy - ip x + p y ) = ^(x + iy - d x - id y ) = -z - 
The only non-vanishing commutators betwen these operators are 

[a,a+] = l , [6,6+] = 1, 

and a + a and 6 + 6 do form a complete set of commuting observables. The wavefunctions ifj 
associated to the lowest Landau level obey aip = 0. They can be labeled by the eigenvalue 
i of the angular momentum L = b + b — a + a which, for the lowest Landau level, reduces to 
6 + 6. Starting from ifj defined by 

atfj = , 6V>o = , 

the normalized solution of which is 



V>0 = -ie 



we can obtain the other states of the lowest Landau level by applying the creation operator 
6 + , with the result 

fa = -^=(b + )% = -^=z e e~^ . (3.2) 



a It is convenient here to use sign conventions different from the ones in ref. 7. 
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For /j, = 0, the Boltzmann factor (3.1) is, up to a multiplicative factor, the squared Slater 
determinant built with the orthogonal wavefunctions ipe,0 < £ < N — 1. 

In the general case /i ^ 0, the Boltzmann factor (3.1) can again be expressed as a 
squared Slater determinant, built with orthogonal wave functions, for the same magnetic 
problem as above. Now the normalized wavefunction -00 has to be replaced by 

\J n cosh n 

Since V'o can be expanded in powers of tanhjU as a linear combination of the previous 
functions (3.2), it is another lowest Landau level wavefunction, which obeys 

aipo = . 

The operators a and a + are unchanged. However, for generating the other mutually 
orthogonal functions which enter the Slater determinant, we must replace b and b + by b 
and b + such that 

= (3.4) 

and such that the commutation relations be preserved. This is easily achieved by a Bo- 
golyubov transformation, i.e. writing for b a linear combination of b and b + , the coefficients 
of which are fixed (up to an irrelevant phase) by the conditions (3.4) and [6, b + ] = 1. One 
finds 

11 

b = (cosh / u)6 + — (sinh(u)6 = (cosh n)(-z — d?) — (sinh / u)(-^ + d z ) (3.5) 

and the orthogonal wave functions are 

4> n = -L(fc+)> . (3.6) 
Vn! 



Using (3.3) and (3.5) in (3.6), one readily sees that the ip n are of the form 

$ n = P n (z)e-^ zS+ ^ tanh » )z2 

where P n (z) is a polynomial of degree n defined by the recurrence relation 

1 r 1 dP (z)i 

P n+1 (z) = -== —^ z P n (z) - (sinh^)^^ (3.7) 
y/n + 1 LcoshjU dz J 
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and Pq{z) = (ncoshfj,) 1 / 2 . A simple change of variable relates (3.7) to the recurrence 
relation 8 obeyed by the Hermite polynomials 6 H n {x), 

H n+1 (x) = 2xH n (x) - 

and Hq(x) = 1, and one obtains 

P n (z) = ^(^tanh^^^TTCosh^)- 1 /^^^^) 

= (7rn!)- 1/2 (cosh^)- n -^" + ... . 

Therefore, by making linear combinations of the columns in the determinant A of (3.1), 
one can replace each z n by (7m!) 1 / 2 (cosh / u) n+ iP n (^), and (3.1) becomes a constant times 
the squared modulus of a Slater determinant built with orthonormal wavef unctions. This 
is another derivation of the result obtained by Di Francesco et al 

4. Surface correlations 

Thus, the density correlation functions associated to the Boltzmann factor (3.1) are 
the same ones as for N fermions occupying the orthonormal orbitals ip n (n = 0, 1, . . . N — l). 
All the n-point densities can be expressed in terms of the projector 

N-l 



^(r,r')=Ji(r)t(r') 



n=0 

= -a(r)a(r') V 2 ~"H n ( - I )H n { - t ) (4.1) 
(cosh^)n! n Vsinh2// "Vsinh2// V ' 

where 

In particular, the density at r is Kn(t, r) and the two-point truncated density (correlation 
function) is 

p T N (r,r') = -\K N (r,r')\ 2 . (4.2) 

b We define the Hermite polynomials as in ref.8, with a weight function exp(— x 2 ), while 
Di Francesco et al used exp(— x 2 /2). 



The particles form an elliptical blob. The correlations are expected to be long-ranged 
only near the boundary 9 ' 10 and we can define a surface density-surface density correlation 
by integrating the long-range part of p%(r, r') with respect to the components of r and r' 
normal to the boundary of the blob. Our aim is to evaluate that correlation function for 
macroscopic distances. 

It is appropriate to work with elliptic coordinates (£, rj) related to the Cartesian co- 
ordinates (x, y) by 

x + iy = [(2N - 1) sinh 2p) 1/2 cosh(£ + irj) . (4.3) 

This choice ensures that the boundary ellipse (see Section 2) is one of the coordinate lines, 
£ = £h, where e~ 2 ^ b = tanhjU, up to a negligible correction of order 1/N (using 2N — 1 
rather than 2N in the definition (4.3) will turn out to be convenient in the following). The 
rescaling factor of the conformal transformation (4.3) is 

HZ, V) = I 1 1] I = m - 1) sinh 2/i] ^ | sinh(£ + irj)\ , (4.4) 
and the surface charge-surface charge correlation is 

poo poo 

(a( V )a( V ')) T = e 2 dthfor,) d?h(?, rf)p% b \r, r') , (4.5) 
Jo Jo 

where p^ b \r, r') is the asymptotic form of pjj(r, r') when |r — r'| is a macroscopic distance 
(it is non- vanishing only if r and r' are near the boundary). 

The surface charge- surface charge correlation is particularly interesting as its asymp- 
totic form is expected to depend on the shape of the plasma, but to be otherwise universal 
(i.e. independent of the microscopic detail). 10 For an ellipse-shaped plasma this form is 11 

! Wv)<rW)) T ~ - 7^rr„ -zJrr, -TT77 ~R • ( 4 - 6 ) 



2nk B T x y " y '" 2TT 2 \e ir i - e ir >'\ 2 h(^ b ,ri)h(^ b ,ri') 

The exact result (4.2) therefore provides an opportunity to test (4.6). 

To evaluate (4.5) and thus test (4.6), our major task is to compute p^ b ^ . For our 
overall strategy in this calculation, we take guidance from the calculation of Choquard 



et al, 12 who calculated p^f b \r,r') in disk geometry. Also, we make use of the known 13 
large- n asymptotic expansion for the Hermite polynomial H n (x) : 

2n/2-l e [n+x 2 -x(x 2 -2n-l) 1 / 2 ]/2 n ] 



H n (x) 



7r l/2 n (n+l)/2 (_^L_ _ ^1/4 



X 



— + (— 1)V2" 

+ 1) 1 / 2 { 2n + l ) J 



.(2n + l)y' 2 ' "2n + l 
We begin by using the integral representation 6 



n+l/2 



H n (x) = -== / (x + zt) n e _t eft 

V7T 



twice in (4.1) to obtain 
^iv(r, r') = 

a(r)a(r / ) 1 



7T COsh 

where 



-= / dt 2 e- t2 — / dt ie -*i V 



JV-1 



{2(tanh//)(Z + it x ){Z' + it 2 )Y 



Z = 



^sinh 2jj 

Following Choquard et al, 12 we make use of the asymptotic approximation 

N-l „ 



u 



u 



N 



to nl (N-l)lu-(N-l) 



i + o( 



N-l 



(u-(N-l)) 



which gives 



-rr / / \ 4a(r)a(r / ).l , . N 
K N (r, r) ~ — v ' v ^ (-tanh^)^ 



1 



>JV-i r°° 



x 



(N-l)l v^F J_ 



x 



7T cosh (I 2 
dt 2 e-^(Z' +it 2 ) 



nN-1 r°° 
N-l z 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



(Z' + it 2 )(Z + ih) 



dtie'** (Z + it!) 



N-l 



(4.12) 



2(tanh//)(Z + + it 2 ) — (AT — 1) 

(as in ref. 12, we can show that the terms ignored in going from (4.9) to (4.12) are 0(1/N) 
relative to the terms kept). 

Consider the double integration in (4.12). Since iV is large, by the method of steepest 
descent, the slowly varying factor on the last line can be replaced by its value at that 



saddle point of each integration which maximizes the modulus of the integrand on the 
second line. Now, the saddle points of the integration over t\ occur when 



f'(t) = 0, where /(£) = t 2 - (N - 1) ln(Z + it) , (4.13) 

which gives 

-Z±(Z 2 -2iV + 2) 1 / 2 
t±(Z) = 2i 

It can be shown that, for Z describing a point near the boundary ellipse, taking t+(Z) 

maximizes the modulus of the integrand (with (Z 2 — 2N + 2) 1 / 2 defined as the positive 

square root when Z is real and larger than (2N — 2) 1 / 2 , and its continuation). Similarly, 

t + (Z') is the saddle point to be chosen for the integral over t 2 . Taking the slowly varying 

factors outside the integral, with ti and t 2 so specified, we see from (4.8) that the remaining 

integrals are precisely the Hermite polynomials, and so 

-r. / is a(r)a(r') .1 , sw 1 
K N (r,r') v ' I ^ (-tanh /j,) N 



yrcosh^ v 2 ^ (N — 1)! 
[Z> + (Z' 2 -2N + 2fl 2 \\Z + (Z 2 -2N + 2) 1 / 2 ]H N _ 1 (Z')H N _ 1 (Z) 
X (| tanh n)[Z' + (Z' 2 — 2N + 2) 1 / 2 ][Z + (Z 2 — 2N + 2) 1 / 2 ] — (N — 1) ' [ ' J 

From (4.2), to compute p N , we seek the large- N behaviour of the squared modulus of 
(4.14), with r and r' in the neighbourhood of the boundary ellipse. Let us consider the 
squared modulus of the product of the Hermite polynomials first. Since, from (4.3) and 
(4.10), 

Z = V2N - 1 cosh(£ + irj), 

(4.7) gives 

2 2AT-6 [(iV _ 1)! ]4 e 2(Ar-l) 1 1 



\Hn-i(Z')Hn-i(Z) 



n 2 (N - 1) 2N | sinh(^ + irj)\ \ sinh(f + irf)\ 

x | e (Ar-|)[e- ( « +i ^cosh(^+^)+C+ir,]|2| e (7V-i)[ e -(«'+ i '' , )cosh^'+ir,')+C'W]|2 _ ( 4 _ 15 ^ 



Furthermore, since 



Z+(Z 2 -2N + 2) 1 ' 2 = (2N - 2) 1 ' 2 e i+l " 1 + O(^) 



the squared modulus of the remaining terms in the second line of (4.14) equals 

4 e 2(£+?') i 

+ 0(-). (4.16) 



|eM-€'+<fa-»»')tanh^- 1| 2 K N' 



We also need 

|a(r)a(r')| 2 = exp [ - (AN - 2)(sinh/j)(e _At cosh 2 £ cos 2 rj + e M sinh 2 £ cos 2 ?]) 



x exp 



(4iV - 2)(sinh/j)(e" M cosh 2 f cos 2 77' + sinh 2 f cos 2 r/)l . (4.17) 



Substituting (4.17), (4.16), and (4.15) as appropriate in the squared modulus of (4.14), 
using Stirling's formula, and substituting the result in (4.2) gives 

T(6)/ / \ tanh/j 1 e 2 ^ +g,) 

Pn [T ' r) ~ Sn 3 (N - 1) cosh 2 fi | sinh(£ + irj) sinh(£' + it/')| \e^'+ i ^') tanh/j - 1| 2 

where 

<7m(£' v) = hitanh/j + 1 — 4(sinh/j)(e _M cosh 2 £ cos 2 77 

+ sinh 2 £ sin 2 77) + e~ 2i cos 2t? + 2f . (4.19) 

Analogous to the behaviour of for the disk, 12 we expect the last exponential in 
(4.18) to be strongly peaked for r and r' in the neighbourhood of the boundary of the 
ellipse, 

i = i h where e" 2 ^=tanh/j. (4.20) 

This is indeed the the maximum of (4.19) occurs when £ = Furthermore, 

9v(£.b, ?]) = is independent of 77. Expansion around this point gives 

Thus, after replacing by in the prefactor of (4.18), we obtain 

T(b) tanh/j e 4 ^ b 1 

Pn ~ ~87r 3 (A^- l)cosh 2 // Isinh^ + ^sinh^ + ^OI |e"> -e^'| 2 



x e -2(t;-t; b ) 2 m b ,v)] 2 e -m' -u'imW)] 2 _ ( 4 _ 2 i) 
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The task of computing p N ( is now completed. To check (4.6), we substitute (4.21) 
in (4.5) and note that since the integrands are sharply peaked at £, £' = we can set 
h(£iV) = = h(£b,v) an d remove these factors from the integrand. Computing the 

remaining Gaussian integrations, using (4.4) (where 2iV— 1 can be approximated by 2N—2) 
and (4.20) to simplify the resulting expression, and recalling that here e 2 /27r/cB^ = 2, we 
see that the predicted form (4.6) is indeed realized. 
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